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We report global analysis results of experimental data for nuclear structure- function ratios j and proton- 
nucleus Drell-Yan cross-section ratios f^'^/cr^y in order to determine optimum parton distribution functions 
(PDFs) in nuclei. An important point of this analysis is to show uncertainties of the distributions by the Hessian 
method. The results indicate that the uncertainties are large for gluon distributions in the whole x region and 
for antiquark distributions at x > 0.2. We provide a code for calculating any nuclear PDFs at given x and 
for general users. They can be used for calculating high-energy nuclear reactions including neutrino-nucleus 
interactions, which are discussed at this workshop. 



1. Introduction 

Although parton distribution functions (PDFs) 
in the nucleon are now determined relatively well 
in the wide range of x, their nuclear modifica- 
tions are not determined accurately. There are 
experimental measurements about nuclear effects 
on the PDFs, for example nuclear data, so 
that their gross properties are known. However, 
the details are not still determined because avail- 
able experimental data exist for a limited number 
of high-energy nuclear reactions. 

On the other hand, accurate nuclear parton dis- 
tribution functions (NPDFs) |ll2IHI4j are needed 
for describing any high-energy nuclear reactions 
such as heavy-ion and neutrino reactions. The 
major purpose of this workshop is to discuss 
nuclear effects on neutrino reactions jSj- Since 
neutrino-oscillation experiments |S] become more 
and more accurate, it is now important to take 
nuclear medium effects into account. 

Current oscillation experiments are done in a 
medium-energy region, so that there are various 
nuclear effects which contribute to the neutrino 
reactions. Among them, we investigate nuclear 
medium effects on the PDFs. Of course, the par- 
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ton distributions are supposed to be used in the 
deep inelastic (DIS) region, so that an extension 
to a smaller region becomes important in or- 
der to use them for the present neutrino reactions. 

In this paper, we report our recent studies on 
the NPDFs by analyzing the nuclear data on the 
structure- function ratios F-^/F^ and Drell-Yan 
cross-section ratios 

'^w/o-w 0- We had ab 
ready reported the NPDF studies on the first ver- 
sion ^ at the NuIntOl workshop and prelimi- 
nary studies after the first version at the Nulnt02 

The most important point in Rcf. 2 is that 
uncertainties of the NPDFs are estimated by the 
Hessian method. For calculating any nuclear re- 
actions, it is especially important to show the un- 
certainties which come from the PDFs. Further- 
more, Drell-Yan and HERMES data are added 
into the data set and the charm-quark distribu- 
tions are included in the new analysis. 

This paper consists of the following. In section 
121 a method is explained for analyzing F2 and 
Drell-Yan data in order to obtain the optimum 
NPDFs. Analysis results are shown in section 13 
We provide a code for calculating the NPDFs at 
given X and Q^, and it is explained in sectional 
The results are summarized in section |5| 
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2 ANALYSIS METHOD 



2. Analysis method 

The parton distribution functions are expressed 
in general by two variables, x and Q^. In lepton- 
nucleus scattering, the is defined by the mo- 
mentum transfer q: = —q^, and the Bjorken 
scaling variable x is given hy x — / {2AIv) with 
the nucleon mass M and the energy transfer v. 
The NPDFs are provided by an analytical form 
at a fixed point Practically, a NPDF 

could be expressed by the corresponding nucle- 
onic distribution multiplied by a function Wi. 



f,^{x,Ql)^w,{x,A,Z)Mx,Ql). 



(1) 



We call Wi a weight function, which indicates nu- 
clear modification. The function is expressed by 
a number of parameters: 
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ai{A, Z) + biX + Cix"^ + diX^ 
^ (l-x)ft • 



(2) 



The parameters are determined by a analy- 
sis with experimental data on structure functions 



and Drell-Yan cross sections a 



pA 
DY- 



Because 



of baryon-number, charge, and momentum con- 
servations, three parameters can be fixed. The 
cubic functional form of the numerator is moti- 
vated by the x dependence of typical data for 
F^/F^, and the factor 1/(1— a:)''* is to reproduce 
the Fermi-motion part at large x. The nuclear 
dependence 1 — l/A^^^ is introduced in Ref. [S] 
simply by considering nuclear volume and surface 
contributions to cross sections. Because different 
physics mechanisms contribute in each x region, 
the overall 1/A^^^ dependence would be too sim- 
ple to describe the nuclear modifications. For the 
NPDFs ff-, we take u^, d^, q^, and by as- 
suming flavor symmetric antiquark distributions 
although they are not symmetric in the nucleon 

nni- 

The parameters are determined by a analy- 
sis with experimental values. The is given by 



/ T}data TDtheo\2 



(3) 



where af*^ is an experimental error, and Rj in- 
j •> 

dicates a ratio, F^ jF^' or cr^'^/o'^^. Leading- 
order expressions are used in the theoretical cal- 
culations. 

The actual calculation is done by running the 
subroutine MINUIT. The subroutine also produces 
a Hessian matrix Fl which has information on pa- 
rameter errors. Using the matrix, we can calcu- 
late the uncertainty of a NPDF: 



\^!\x)f^^x^Y. 



{ dj\x,i) 



X Rl 



df^{x,0 



(4) 



where is a parameter, ^ indicates the optimum 
set of the parameters, and 5f^{x) is the uncer- 
tainty of the nuclear PDF f^{x). The Ax^ value 
determines a confidence region and it depends on 
the number of parameters. 

The kinematical region of the experimental 
data for F^ and a^j^ is shown in Fig. ^ In com- 
parison with the nucleon data, the x range is still 
limited. Namely, the small-x {x = 10~^ — 10"'^) 
data are not taken. The Drell-Yan data are taken 
in the large region. There are 606 data points 
for the F^/F^ type, 293 for F^/F^ [A' ^ D), 
and 52 for the Drell-Yan. The total number of 
data is 951. These data are taken for the targets: 
deuteron (D), helium-4 (^He), lithium (Li), beryl- 
lium (Be), carbon (C), nitrogen (N), aluminum 
(Al), calcium (Ca), iron (Fe), copper (Cu), kryp- 
ton (Kr), silver (Ag), tin (Sn), xenon (Xe), tung- 
sten (W), gold (An), and lead (Pb). 
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Figure 1. Kinematical region is shown for exper- 
imental data. 



3.1 Comparison with experimental data 
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3. Results 

3.1. Comparison with experimental data 

The NPDFs in Eqs. ((TJ and Q are given at 
(5q=1 GeV^, and they are evolved to the exper- 
imental points for calculating the total in 
Eq. 0. The optimum parameters are obtained 
by minimizing x^- The minimum becomes 
1489.8 for the 951 data. The uncertainties of the 
NPDFs are estimated with Ax^=10.427 in order 
to show the one-cr-error range. It is chosen by 
considering that the number of the parameters is 
nine piTT] . 

Each contribution is listed in Tabled There 
is a tendency that the x^ values are large for small 
nuclei. The Li/D, Be/D, and C/D ratios have 
X^ values, 88.7, 44.1, and 130.8, for the number 
of data points, 17, 17, and 43 points, respectively. 

Table 1 



Each contribution. 



nucleus 


# of data 


x' 


"He/D 


35 


56.0 


Li/D 


17 


88.7 


Be/D 


17 


44.1 


C/D 


43 


130.8 


N/D 


162 


136.9 


Al/D 


35 


43.1 


Ca/D 


33 


42.0 


Fe/D 


57 


95.7 


Cu/D 


19 


11.8 


Kr/D 


144 


126.9 


Ag/D 


7 


12.8 


Sn/D 


8 


14.6 


Xe/D 


5 


2.0 


Au/D 


19 


61.6 


Pb/D 


5 


5.6 


F^/F^ total 


606 


872.8 


Bc/C 


15 


16.1 


Al/C 


15 


6.1 


Ca/C 


39 


36.5 


Fe/C 


15 


10.3 


Sn/C 


146 


257.3 


Pb/C 


15 


25.3 


C/Li 


24 


78.1 


Ca/Li 


24 


107.7 


F2^/F2^ total 


293 


537.4 


C/D 


9 


9.8 


Ca/D 


9 


7.2 


Fe/D 


9 


8.1 


W/D 


9 


18.3 


Fe/Be 


8 


6.5 


W/Be 


8 


29.6 


Drell-Yan total 


52 


79.6 


total 


951 


1489.8 



Medium and large size nuclei are generally well re- 
produced except for the Sn/C and Ca/Li ratios. 
The Drell-Yan data are also well explained except 
for the W/Be ratios. 

Typical results are shown in Figs. |3and|21 In 
Fig. El the ¥2°" /F2 data are compared with the 
fit result at (5^=5 GeV^. The shaded area is the 
uncertainty range due to the NPDF uncertainties 
at Q^=5 GeV^. One should note that the experi- 
mental data are taken at various points which 
are not equal to 5 GeV^. Therefore, the curve 
cannot be, strictly speaking, compared with the 
data. However, considering that the depen- 
dence is small in general, we find that the data 
are reproduced well by the parametrization. 

The Drell-Yan data (^^^y I'^^dy compared 
with the fit result in Fig. 13 The parametrization 
curve and the uncertainty range are calculated at 
(5^=50 GeV^, whereas the data are taken at var- 
ious points. In the region x < 0.1, the cross- 
section ratio is almost the same as the antiquark 
ratio <f°'/q^, so that the data play a role of fixing 
the nuclear antiquark distributions at a; ^ 0.1. 

The actual comparison with the experimental 
data should be done at the same points. In 
order to estimate the fit result, we show the ratios 
^^exp _ j^theoy j^theo ^ig. H Here, indi- 
cates an experimental F2 ratio and R*'^^° does a 
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Figure 2. Comparison with the F!^"" /F2 data. 
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Figure 3. Comparison with the <J^y /'^dy data. 
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3 RESULTS 
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Figure 4. Fit results are compared with exper- 
imental data for R = F^/F^ at the same ex- 
perimental points. The fractional differences 

(^exp _ j^theoy j^theo g^^^ ghowu. 



ratio by the parametrization. The theoretical ra- 
tios are calculated at the experimental points. 
Among the used data, only F^ / F^ type data are 
shown in Fig. 0] In general, the data are well ex- 
plained by the parametrization. However, there 
are slight deviations in small nuclei as indicated 
in Tabled For example, the lithium data at small 
X are underestimated. On the other hand, the tin 
data are overestimated at small x. We need more 
complicated A dependence for explaining all the 
nuclei. 

3.2. dependence 

Next, Q"^ dependence of F2 is calculated and it 
is shown with the F^ /F.p data in Fig. |S1 where 
A'' indicates nitrogen. The figure shows that the 

dependence of the ratio F^ /F^ is not very 
obvious experimentally, which makes the determi- 
nation of the nuclear gluon distributions difficult. 
The ratio tends to decrease with increasing at 
X = 0.035 and x = 0.045. There is a same ten- 
dency in the F^''/F^ data by the HERMES col- 
laboration. However, the ratio i^2^"/i^2^ increases 
with increasing at a; = 0.0125 ~ 0.045 accord- 
ing to the NMC collaboration [2]. It seems that 
the variations are not consistent each other be- 
tween the HERMES and NMC data. Since such 

dependence is crucial in determining nuclear 
gluon distributions, we hope that the depen- 
dence will be accurately measured in the small- a; 
region. 









.Y=0.035 




Figure 5. dependence of F^ /F^ is shown in 
comparison with the data. The curves indicate 
fit results. 
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3.3. NPDFs with uncertainties 

From the fit, we obtain the optimum NPDFs 
with the uncertainties. As an example, the weight 
functions are shown for the calcium nucleus in 
Fig. These functions indicate nuclear mod- 
ifications of the valence-quark, antiquark, and 
gluon distributions in the calcium at Q^=\ GeV^ 
by definition. The shaded areas indicate the un- 
certainties estimated by the Hessian method. 

The valence-quark distributions are well de- 
termined at medium x by the F2 data because 
the is dominated by them in this x region. 
In the small-x region, the valence distributions 
are fixed by baryon-number and charge conserva- 
tions. These are the reasons why the uncertain- 
ties are rather small. However, it is worth while 
investigating the valence-quark shadowing at the 
NuMI and neutrino-factory projects ^2] despite 
the present result. 

The antiquark distributions are determined at 
small X because of the shadowing data of F2 at 
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Figure 6. Nuclear modifications of valence-quark, 
antiquark, and gluon distributions in the calcium 
nucleus at (5^ = 1 GeV^. 



small X. They are also fixed at a; ~ 0.1 by the 
Drell-Yan data. We notice that the uncertainties 
are fairly small in these regions. However, they 
cannot be determined at x > 0.2 as indicated by 
the large uncertainties in Fig. We need new 
experiments such as J-PARC and Fermilab-P906 
in this large- a; region 

The gluon distributions cannot be fixed at this 
stage although they seem to be shadowed at small 
X. In fact, the uncertainties are huge in the whole 
X region. The difficulty is mainly because accu- 
rate scaling violation data are not available at 
small X as shown in Fig. [3 The variations 
in the x ^ 0.01 region are not measured accu- 
rately, and it is also obvious from Fig. ^ that 
the small- a; (x = 10"'^ ^ 10~^) data themselves 
do not exist. In addition, a next-to-leading-order 
analysis could be important for incorporating the 
gluon contributions into the parametrization. 

From the analysis, we obtain the NPDFs from 
the deuteron to a heavy nucleus with A ^ 208. 
We provide the NPDFs for general users by pro- 
viding a code, which is explained in the next sec- 
tion. The NPDFs could be used for any high- 
energy nuclear reactions. In order to use them 
for the present neutrino oscillation experiments 
in a medium-energy region, an extension from the 
DIS region to the resonance one should be investi- 
gated. Such an effort has been done, for example, 
in Ref. for the structure function F2 of the 
nucleon. Obviously, we need a similar study for 
extending the NPDFs to the smaller region in 
order to calculate nuclear corrections to the cross 
sections in the neutrino oscillation experiments. 

4. Code for calculating NPDFs 

A useful code could be obtained from the web 
site in Ref. P'. In the package (npdf04.tar.gz), 
grid data, npdf04 subroutine (npdf04.f), and a 
sample file (sample. f) are provided. Running the 
subroutine, one obtains the NPDFs at given x 
and for the analyzed nuclei, D, "'He, Li, Be, 
C, N, Al, Ca, Fe, Cu, Kr, Ag, Sn, Xe, W, Au, 
and Pb. The kinematical ranges for the provided 
NPDFs are 10"^ < a; < 1 and 1 GeV^ < < 
lO'* GeV^. The smallest x point of the data is 
a^miri^O.OOSS, so that one should note that the 
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5 SUMMARY 



NPDFs at X < x„iin are not tested experimen- 
tally. However, we provide the small-x distribu- 
tions in case that one uses the NPDFs for inte- 
grating them over a wide range of x. The maxi- 
mum Q2 of the data is 173 GeV^ The vari- 
ations above this point are also not tested ex- 
perimentally. We simply assumed the standard 
DGLAP evolution equations for extending the 
region up to = 10* GeV^ for feasibility 
studies of future experimental facilities. 

For a nucleus other than the provided ones (D, 
He, • ■ •, Pb), the NPDFs could be calculated by 
following the procedure in Appendix of Ref. |2] . 
Namely, the nuclear dependent parameters a^^ , 
ad^, ag are calculated by Eq. (Al) of Ref. and 
other parameter values are taken from Table H. 
Next, the NPDFs for a nucleus can be calculated 
at (5^=1 GeV^ by using the analytical expressions 
with the determined parameters. Then, one may 
evolve the NPDFs to a given point by one's 
evolution code or, if it is not available, by the 
code in Ref. I15j . The analyzed nuclei are up 
to A ~ 208. However, variations from the lead 
NPDFs to those of the nuclear matter {A oo) 
are small, so that one may use the analysis results 
for large nuclei with A > 208. 

5. Summary 

We have investigated a parametrization of 
NPDFs by analyzing experimental data on 
and Drell-Yan processes. In addition to obtaining 
the optimum NPDFs, we calculated their uncer- 
tainties by the Hessian method. The antiquark 
distributions at small x and valence-quark distri- 
butions are well determined by the data. How- 
ever, the antiquark distributions at x > 0.2 and 
the gluon distributions cannot be determined re- 
liably. We obviously need future experimental 
efforts for determining these distributions. One 
could use the obtained NPDFs by running the 
provided code. 
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